Two-photon and three-photon blockades in driven nonlinear systems 
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Photon blockade, in analogy to Coulomb's or phonon blockades, is a phenomenon when a single photon in a 
nonlinear cavity blocks the transmission of a second photon. This effect can occur in Kerr-type systems driven 
by a laser due to strong nonlinear photon-photon interactions. We predict the occurrence of higher-order photon 
blockades where the transmission of more than two photons is effectively blocked by single- and two-photon 
states. This photon blockade can be achieved by tuning the frequency of the laser driving field to be equal 
to the sum of the Kerr nonlinearity and the cavity resonance frequency. We refer to this phenomenon as two- 
photon blockade or two-photon state truncation via nonlinear scissors, and can also be interpreted as photon- 
induced tunnelling. We also show that, for a driving-field frequency fulfilling another resonance condition and 
for higher strengths of the driving field, even a three-photon blockade can occur but less clearly than in the 
case of single- and two-photon blockades. We demonstrate how various photon blockades can be identified 
by analyzing photon-number correlations, coherence and entropic properties, Wigner functions, and spectra of 
squeezing. We show that two- and three-photon blockades can, in principle, be observed in various cavity and 
circuit quantum electrodynamical systems for which the standard single-photon blockade was observed without 
the need of using two-photon driving fields or Kerr media exhibiting higher-order nonlinear susceptibility. 

PACS numbers: 42.50.Dv, 42.50.Gy, 42.50.Lc 



I. INTRODUCTION 

In nonlinear optical systems driven by a coherent classical 
light in a cavity, a single photon can impede the transmission 
of other photons. This phenomenon is referred to as photon 
blockade (PB) [1], in close analogy to the phenomenon of 
Coulomb's blockade |2], where electron transport is blocked 
by a strong Coulomb interaction in a confined structure. In 
PB the next photon can enter the cavity only if the first photon 
has left it; thus, a sequence of single photons can be generated 
and such system can act as a single-photon turnstile device. 
Evidently, the PB changes classical light into highly nonclas- 
sical light exhibiting, in particular, photon antibunching and 
sub-Poisson photon-number statistics. 

The PB can be interpreted as nonlinear optical-state trunca- 
tion or nonlinear quantum scissors , since the infinitely- 
dimensional Fock-state expansion of a classical driving field 
is truncated at the single-photon Fock's state. The required 
nonlinearity can be induced by a strong interaction between 
cavity and two-level (natural or artificial) atom. 

The PB has been predicted in various setups in cavity quan- 
tum electrodynamics (QED) jsUTTIl. and recently also in cir- 
cuit QED dill. The PB was first demonstrated experimen- 
tally with a single atom trapped in an optical cavity vi\M . This 
experiment was considered "a landmark event in the field of 
quantum optics and laser science" llSil . In solid-state systems, 
the PB was experimentally demonstrated with a quantum dot 
in a photonic crystal cavity (l^ and with a single supercon- 
ducting artificial atom coupled to a microwave transmission- 
line resonator (superconducting "cavity") lll2[ Hsll . The PB 
was also predicted in quantum optomechanical systems ITtIi . 



An analogous phenomenon of phonon blockade was predicted 
for an artificial superconducting atom coupled to a nanome- 
chanical resonator ifisll . The phonon blockade can be detected 
via the PB if this system is additionally coupled with, e.g., a 
superconducting microwave cavity liioll . 

The PB can have applications in quantum state engineer- 
ing for the controllable generation of a train of single pho- 
tons exhibiting highly nonclassical photon statistics. This sug- 
gests that the PB (together with phonon blockades) can also 
be used as an indicator of nonclassicality of mechanical sys- 
tems iflsllTgll . As another example of quantum engineering, 
the PB was also studied in the context of cavity electromag- 
netically induced transparency in, e.g., Refs. [6, 9]. The pos- 
sibility of tunable (by a classical driving field) transmission 
from photon blockade to photon transparency in circuit-QED 
systems was described in Ref. f^^l- 

Moreover, single-photon tunnelling, in close analogy to 
single-electron tunnelling, was experimentally observed dur- 
ing light transmission through individual subwavelength pin- 
holes il Id . This effect was explained in terms of the PB, 
analogously to the Coulomb blockade demonstrated in single- 
electron tunnelling. A close analogy between the PB and 
Coulomb blockade was also analyzed in Ref. |20] by present- 
ing, e.g., a photonic analog of the Coulomb staircase. Such 
studies shed more light on quantum simulations of condensed- 
matter phenomena via optical effects and vice versa. 

In this paper, we show the possibility of observing of multi- 
photon blockade in the standard cavity QED and circuit QED 
systems, where the single-photon blockade was already ob- 
served. These systems consist of a Kerr nonlinearity in a cav- 
ity driven by a classical weak field. We will show how to 
change the intensity and frequency of the driving field for a 
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FIG. 1. (color online) Schematic energy-level diagram explaining 
the occurrence of the fc-photon blockade (and the fc-photon-induced 
tunneling) in terms of the fc-photon transitions induced by the driving 
field satisfying the resonance condition — 0, which corresponds 
to the driving-field frequency uja = i^;^''' — luq + x{k — 1). Due to 
the Kerr-type nonlinearity (induced by, e.g., a qubit), the cavity-mode 
levels isi"' = nhujo (shown on the left) become non-equidistant as 
A_E„+i,„ = En+i — En / const, where E„ fa n[&Jo + (n - l)x] 
(for n — 0,1, ...) are the eigenvalues of the Hamiltonian H, given 
by Eq. (TJ, assuming e <^ x- 



given cavity frequency and the Kerr nonlinearity in order to 
observe the two- and three-PBs. 

Those multi-PBs can also be explained as photon-induced 
tunnelling. That is, when there is a photon inside the cavity, 
then the second or third photon can be absorbed by the cavity 
via two-photon or three-photon processes. 

Effects leading to multi-PB were already studied in the lit- 
erature as generalizations of the single-PB. For example, sev- 
eral systems based on higher-order parametric driving pro- 
cesses and/or higher-order Kerr nonlinearity were analyzed in 
Refs. ll2ll - l23|] for the observation of the nonstationary-field 
multi-PB. In the Conclusions, we will discuss shortly the for- 
mal differences and crucial experimental advantages of our 
approach over these methods. 

There is another approach to the PB, which is based on lin- 
ear systems, while the nonlinearity is induced by measure- 
ments. This method is usually referred to as the linear opti- 
cal state truncation or linear scissors fl^. The multiphoton 
state truncations (or multi-PBs) implemented via linear scis- 
sors were studied in Refs. |25, 26], as a generalization of the 
single-photon state truncation [24,, i27il . It is worth stressing 
that the PBs studied in this article are based on completely 
different principles and resources. We use, as in the original 
single-PB proposals fH [H, nonlinear systems (which cause 
nonlinear photon-photon interactions) without measurements. 
Thus, the method used here can be referred to as nonlinear 
scissors [4]. 

The paper is organized as follows. In Sec. II, we explain 
the occurrence of the single-PB and describe the method to 
observe blockades up to two and three photons. In Sec. Ill, 
we demonstrate analytically the occurrence of the two-PB in 
comparison to the single-PB. In Sec. IV, we discuss vari- 
ous signatures of the two-PB revealed by the photon-number 
statistics, entropies, Wigner functions and spectra of squeez- 
ing. We summarize our main results in the concluding section. 
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FIG. 2. (color online). Photon-number probabilities P„ — {n\pss\n) 
for the steady-state solutions pss of the master equation for the 
Hamiltonian H^l!;^ (0) as a function of the tuning parameter k assum- 
ing the driving strengths (a) e = 67 and (b) e — 11. 567. Moreover, 
we assume: the Kerr nonlinearity x ~ 30"f, the damping constant 
7 = 1, and the mean number of thermal photons nth ~ 0.01. Fig- 
ures also show the truncation fidelity Fm (for m — 2,3) and the Fano 
factor F. Note that the field exhibits sub-Poisson (super-Poisson) 
photon-number statistics if F < 1 (F > 1). It is seen in both figures 
that resonances at fc = 1 and 2 can be interpreted as the single- and 
two-photon blockades, respectively. But the three-photon blockade 
at fc = 3 is apparent in (b) only. 



II. FROM SINGLE-PHOTON TO MULTIPHOTON 
BLOCKADES 



A. Hamiltonians 



Photon-blockade effects can be observed in a cavity with 
a nonlinear medium coherently driven by a laser field as de- 
scribed by the following effective Hamiltonian ^ Isj] : 
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FIG. 3. (color online). Photon-number probabilities Pn = {n\pss\n) 
and fidelities Fm = X]™=o ^" m-photon truncation for the 

Hamiltonian -ff^'ot (0) for the resonances at (a) fc = 2 and (b) = 3 
as a function of the driving strength e assuming the Kerr nonlinearity 
X ~ 3O7, the damping constant 7 = 1, and the mean thermal-photon 
number fith ~ 0.01. Note that Po Pi ?a P2 for the driving 
strength e = 11.56 and fc = 3. This value of e is chosen in Figs. 2(b), 
5, 8 (for k — 3), and 11. Also note that in (a) for k = 2, there are 
two crossings, one for Po = Pi and another one for Po = P2. 



where Ai = ~ "^d- The single-PB can be observed in the 
resonant case lu^ — luq assuming that the driving strength e 
is much smaller than the Kerr nonlinearity x- This effect can 
be interpreted as the single-photon Fock state blockading the 
generation of two or more photons. 

We observe that Hamiltonian ^ can be rewritten as fol- 
lows: 

Hlll{l^k) = hAkCi'^a + hxa^ a{a'' a - k) + he{a + 0)), (5) 
where the frequency mismatch is 

Afc = Wo + x{k -I) ~ujd (6) 

for some positive k. For convenience, we shall refer to fc as a 
tuning parameter In the special case of A: = 1, this Hamilto- 
nian reduces to the standard form, given by Eq. (3). 



where d [a)) denotes the annihilation (creation) operator, 
is the resonance frequency of the cavity, ujd is the driving laser 
frequency, x > is the Kerr nonlinearity, i.e., the photon- 
photon interaction strength proportional to the real part of the 
third-order nonlinear susceptibility Re(x''^^), and e > is the 
driving strength (the Rabi frequency of the laser). 

Equation ([T]i presents an effective Hamiltonian, which can 
be obtained from various microscopic Hamiltonians describ- 
ing a variety of systems. For example, one can analyze a 
quantum two-level system (qubit) off-resonantly coupled to 
a driven cavity. This system can be described in the rotating- 
wave approximation by the following Hamiltonian: 

H — —hujaz + hw'^a) a + hg{cT^ a + a) a^) (2) 

where the first three terms correspond to the standard Jaynes- 
Cummings Hamiltonian and the last term, as in Eq. ([TJ, 
describes the interaction between the quantum cavity mode 
and the classical driving field with strength e. Moreover, 
'^z — |e)(e| — \g){g\ denotes the Pauli operator, (t+ — \e){g\ 
(cr~ = Ig) (e|) is the qubitraising (lowering) operator, \g) (|e)) 
is the ground (excited) state of the qubit, uj is the qubit transi- 
tion frequency, and Wq is the resonance frequency of the cav- 
ity in Eq. (|3]l. Moreover, other quantities are explained below 
Eq. ([TJ. This Hamiltonian in the dispersive approximation, 
where the qubit remains in its ground state, can be reduced to 
the Hamiltonian, given by Eq. ^ (for a derivation see, e.g., 
Refs. 1 18, 28] in the circuit QED context). 

The unitary operation [/ = exp(— iwdO^at) transforms the 
Hamiltonian ([TJ into 

H.ot^U^HU -ihlJ^^U. (3) 
at 

Thus, in the rotating-frame, one obtains the following time- 
independent Hamiltonian, 

h[^^ = ftAia^a + hx{a)fa^ + ht{a + a)), (4) 



B. Steady states as a function of the tuning parameter k 

The evolution of the system, given by Eq. ((5j, for the re- 
duced density operator pit) under Markov's approximation, 
can be governed by the standard master equation I29J . which 
for the steady state p^s — p{t — > 00) is given by 

~ Pssaa^) + ^("'th + l)(2a/5ssa^ - a^apss - Pssa^a),(l) 

where 7 denotes the damping constant, nth = 
{exp[?za;/(fcBr)] — 1}^^ is the mean number of ther- 
mal photons, ks is the Boltzmann constant, and T is the 
reservoir temperature at thermal equilibrium. The master 
equation can be given in terms of a Liouvillian superoper- 
ator and the steady-state solution p{t) can be obtained by 
applying, e.g., the inverse power method ll30ll . 
The single-PB can occur if the conditions 

7 < e < X (8) 

are satisfied. Hereafter, we only analyze the resonant case 
Afc = 0, which is related to the resonant fc-photon transitions 
shown in Fig. 1 . This condition implies that the tuning param- 
eter k is related to the Kerr nonlinearity, the driving-field and 
cavity frequencies as follows 

k = {ujd-uJa)/x + l- (9) 

Figure 2 shows how the photon-number probabilities P„ of 
the steady states depend on the tuning parameter k. By analyz- 
ing this figure, one can discover various kinds of PB effects, 
which appear not only for the standard resonant case of fc = 1 
but also for fc 7^ 1. For example, in the special case of fc = 2, 
a higher-order effect occurs with at most two photons effec- 
tively generated in the system. We refer to this effect as the 
two-photon blockade, which means that the single- and two- 
photon Fock states blockade the generation of more photons. 

Figure 2(b), which is similar to Fig. 2(a) but obtained as- 
suming larger driving strength e, clearly shows a three-photon 
blockade for fc = 3, which refers to a phenomenon where the 
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Fock states |m) (for m = 1,2,3) blockade the transmission of 
the Fock states \n) with higher photon number (i.e., n > 3). 
Note that, in Fig. 2 for fc = 3, the probabiHties of generat- 
ing 1 2) and |3) are nonzero but are much lower than the 
probability of observing the vacuum state. So, this behavior 
for the parameters of Fig. 2(a), contrary to Fig. 2(b), cannot 
be considered a genuine three-PB. 

It is also quite evident why it is necessary to increase the 
strength of the driving field to obtain a good-quality three- 
PB since we need a larger average photon number For weak 
driving, the average photon number of the driving field is less 
than three photons, thus this is not enough to induce a three- 
photon transition as shown in Fig. 1. 

In Fig. 2, in particular, we show the fidelity of the m-photon 
truncation defined as 



n\pss\n) 



(10) 



We refer to the to-PB if the truncation fidelity F„i w 1 and 
Fn ^ 1 for n < m. In Fig. 2(a), it is seen that Fi ss 1 at fc = 
1 (^2 « 1 at /c = 2) corresponding to the single-PB (two-PB). 
But the resonance at fc = 3, for the chosen driving strength, 
can hardly be considered the three-PB for the driving strength 
e = 57. By contrast, F3 « 1 at fc = 3 for the much higher 
driving strength (e — 11. 567), as shown in Fig. 2(b). Thus, 
we interpret the latter case as the true three-PB, as already 
mentioned. 

By analyzing the dependence of the photon-number proba- 
bilities Pn ~ {n\pss\n) on the tuning parameter k, as shown in 
Fig. 2(a), we observe that there is a clear dip in the vacuum- 
state probability Pq and a peak in the single-photon proba- 
bility Pi at fc = 1 and 2. The two-photon probability P2 is 
nonzero only near k = 2. The peak maximum in P2 does not 
occur at fc = 1 but it is evident at fc = 2, which implies that the 
single-photon truncation fidelity Fi differs from 1 near k = 2. 
By contrast, the two-photon truncation fidelity F2 on the scale 
of Fig. 2(a) is practically equal to one. Thus, the contribution 
of Fock's states with three or more photons is negligible. For 
these reasons, we refer to the two-PB at fc = 2. One can see in 
Fig. 2(a) a very slight dip in Pq at k ~ 3, which can be, how- 
ever, much deeper for the approximately twice-larger driving 
strength e, as shown in Fig. 2(b). The dip in Pq at k = 3 is 
accompanied by the clear appearance of the peaks of P„ for 
the photon numbers n — 1,2,3. It is seen that for fc = 3 (but 
also close to this point), the three-photon truncation fidelity 
F3 « 1 contrary to F2 ^ 1. This explains why we refer to 
this effect as the three-PB at fc = 3 for a suitably large driving 
strength e as, e.g., in Fig. 2(b). 

Figure 3 shows how by increasing the driving strength e one 
can change the steady-state probabilities of the generation of 
n-photons for chosen values of the Kerr nonlinearity, damping 
constant and for the resonance conditions (a) fc = 2 and (b) 
fc = 3. We also depicted the two- and three-photon truncation 
fidelities, which practically equal to 1 for the all values of e in 
figures (a) and (b), respectively. It is seen that for the small 
driving strength one cannot generate photons in the cavity in 
the steady-state limit, as Pq ~ 1- This is a trivial case, which 
can be interpreted as the zero-PB. For larger values of e, the 



true single- and two-PB effects are observed. In particular, it 
is seen in Fig. 3(b) that for the driving strength e = 11.56 the 
probabilities P„ of observing n = 0,1,2 photons are approx- 
imately the same. The three-photon probability is smaller but 
still nonzero. We have chosen this e to analyze various prop- 
erties of the three-PB in Fig. 2 and other figures. 



III. ANALYTICAL DESCRIPTION OF PHOTON 
BLOCKADES 

A. Steady-state photon blockade 

Here, we explain analytically the standard (i.e., steady- 
state) PB effects, i.e., by including dissipation in the system 
in the infinite-time limit. 

We start from the infinite-dimensional Hamiltonian, given 
by Eq. (|5]l in the resonant case ~ for a given fc, and for- 
mally truncate it to a finite-dimensional Hilbert space. For ex- 
ample, to show explicitly the two-PB (fc = 2), one should an- 
alyze the Hamiltonian at least in the four-dimensional Hilbert 
space: 
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(11) 



which is given in the standard Fock basis. We have to show 

(2) 

that the contribution of three-photon terms {3\pss \n) (forn — 
0, 1, 2, 3) in the steady-state density matrix is negligible. 

The conditions given by Eq. (O should be fulfilled to ob- 
serve a PB effect, so let us denote 
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(12) 



where (5 <C 1 and dS ^ 1. To find an analytical approximate 

(2) 

steady-state solution for pss \ we substitute the truncated four- 
dimensional Hamiltonian into the master equation, given by 

Eq. O, and solve the set of equations for each of the Fock- 

(2) 

state elements of the density matrix (m|pss for m,n = 
0, 1, 2, 3, 4. Then we expand these solutions in power series 
of 5 and neglect terms proportional to 5"^ or higher powers. 
Thus, we find the following approximate steady-state solution 
up to 5^: 



^(2) 



1 



l + 8d2 



/I 



V 



+ 2d2 

x5 

idV2 
z5 



x*S 

yS 




idV2 

y*S 

2d^ 



z*5 \ 


l^/id^i 
IVid'^S 

(13) 

72^2 {2d + i) , and z = 



where x = ~2d^ — 2id'^ + d, y 
^i^/6d^. Our numerical solutions slightly differ from Eq. (ITJt 

for the parameters chosen in Fig. 2(a), but it is clearly seen 

(2) 

that the contribution of the three-photon terms (3|/5ss |3) can 

(2) 

be neglected, contrary to the terms (mjpss Iti) with smaller 
number m of photons. Also other elements {3\pss \m) are 
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either equal to zero (for m = 1) or proportional to 5 <C 1 
(for m — 0, 2), so they can be neglected. Thus, our solutions 
explain the two-PB. 

For comparison, let us analyze the single-PB described by 
the infinite-dimensional Hamiltonians, given by Eqs. (HI or Q 
for fc = 1, but truncated to the three-dimensional subspace, as 
given by 











(1) 



,{0) = I e V2e 
V2e 2x 



(14) 



By performing calculations analogous to the former case, we 
find the following steady-state solution up the S^: 



y6^ -\</2d8 j 



x*S y*5^ \ 



(15) 

where x = {2d + i) and y — |\/2d {d + i), which re- 
duces to 




-{d+U)5 



-{d-\i)5 



\dV2S (16) 
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0.2 0.3 0.4 
rescaled time 5t 

FIG. 4. (color online). Two-photon blockade without dissipation: 
Evolution of the photon-number probabilities P„ for the tuning pa- 
rameter k = 2 according to our precise numerical calculations 
j^pcxact^ thin curves) in the 100-dimensional Hilbert space and ap- 
proximate solutions (P,''^""'^, thick curves), given by Eqs. J17t -( l22t . 
obtained in the truncated four-dimensional Hilbert space. Excel- 
lent agreement of these solutions implies that the truncation fidelity 
F2 ~ Po + Pi + P2 is almost exactly equal to one (as shown by 
black line). This explains the meaning of the optical-state truncation 
or the nonstationary-state two-photon blockade. It is seen that the 
main contribution to F2 is from the Fock states |0) and |2), which 
corresponds to the two-photon transitions shown in Fig. 1 . 



assuming (5 ss 0. It is seen in Eq. (HSll that (2|p£'| n) vanishes 
for n = 0, 2 and is (5 ^ 1 for n = 1. So, the contribution 
of all the two-photon states can be neglected, which clearly 
explains the physical meaning of the single-PB for k = 1. 
Note that (2|pss^ |2) 7^ as described by Eq. (O, which cor- 
responds to the two-PB for k = 2. 



B. Photon blockade without dissipation 

Here, we shortly describe nonstationary-state PB assuming 
no dissipation. The main results of this section are summa- 
rized in Figs. 4 and 5. 

To describe the two-PB, let us formally confine the Hilbert 
space of our system to four dimensions. Thus, we use the 
Hamiltonian, given by Eq. ( fTTT i. Its exact eigenvalues and 
eigenvectors can be calculated analytically, but they are too 
lengthy to be presented here. Instead of this, we define a small 
parameter (5 ^ 1 as the ratio the driving strength e and the 
Kerr nonlinearity x as given in Eq. (fT2] | but, for simplicity, we 
assume here d = 1. Then, we find the power-series expan- 
sion in 5 and keep the terms up to 5^ only. Thus, we find the 
following eigenvalues: 

Ai«-x(3(52 + l), A2~x<5'a;_, 
As^X'^V, A4«x(<52 + 3), (17) 




0.4 0.8 1.2 



rescaled time 5t 

FIG. 5. (color online). Three-photon blockade without dissipation: 
Evolution of the photon-number probabilities Pn as in Fig. 4 but for 
fc = 3 and all the parameters (except 7 = 0) as in Fig. 2(b). It 
is seen that the three-photon truncation fidelity F3 slightly deviates 
from one. For clarity, we present here only precise solutions obtained 
in a large-dimensional Hilbert space. It is evident that solely the Fock 
states |0) and |3) are interchangeably highly populated, which cor- 
responds to the three-photon transitions shown in Fig. I. The proba- 
bilities Pn for n = 1, 2 are relatively small but not negligible, while 
Pn for n > 3 can practically be ignored. 
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and the corresponding eigenvectors: 

|Ai) « iVi[2%/2(5|0) - 2\/2 (3(5^ + 
+4<5|2)- V3<52|3)], 
iV_[((5^ +3)|0) +3a;_(5|l) 
+ {x_d^ ~3)\2) + V3S\3)], 
N+[{d^ + 3)\0) + 3x+5\l) 
-ixj'~3)\2)-V3S\S)], 



1 ID 



IA2 



I A, 



IA4) « Ni[S^\l) + 2V2S\2) + 2\/6|3)], 



(18) 

(19) 

(20) 
(21) 



where x± — 1 ± V2, and the normalization constants are 
iVf ^ « 8(1 + 9^2), N^^ « 6(3 + 5 ± 2V2S^), and iV4-2 sa 
8(3 + S'^). So, assuming no damping and no photons initially 
in the cavity, the time evolution of the truncated system can 
be given by 



\m) 



51 ^^p( 



-^X,t){X,\0)\X,)+O{5^). 



(22) 



Note that, since the initial state is |0), there is no contribution 
of IA4). After substituting Eqs. (fTTt-dSTTi into Eq. ( |22] |. we 
find that 



53 



^(2 - 3e** + m) + 0{S^) = 0{6^), (23) 
v3 

which shows explicitly the negligible contribution from the 
three-photon Fock state in the generated state \ijj{t)). Analo- 
gously we can show no important contribution from the Fock 
states with more photons. Thus, we conclude the occurrence 
of two-PB in the dissipation-free system. 

In Fig. 4, we compare (i) the analytical approximate solu- 
tion, given by Eq. (l22l l. with numerical solutions obtained in 
the Hilbert spaces of dimension (ii) A'dim — 100 (which effec- 
tively corresponds to A'dim — 00) and (iii) Adim = 4 (without 
applying expansions in the power series of S) for some param- 
eters satisfying the conditions On the scale of Fig. 4, there 
is apparently no difference between the solutions (ii) and (iii) 
at all, and very tiny discrepancy between them and the analyt- 
ical approximate solution (i) for the relatively large 5 — 1/6. 
This excellent agreement between the approximate and pre- 
cise solutions, convincingly demonstrate the blockade up to 
the two-photon state |2). It is worth stressing that, although 
the contribution of the three-photon state |3) is very small, the 
calculations have to be performed in the Hilbert space includ- 
ing the state |3). 

For comparison, we recall that the truncated three- 
dimensional Hamiltonian H^^l^-^^(0) under standard condi- 
tions, given by Eq. (O, leads to the following simple evolu- 
tion 

\^{t)) = cos(ei)|0) - isin(et)|l) +0{S^), (24) 

assuming \ip{0)) = |0). As shown by Leoriski and Tanas in 
Ref. [3J, this solution well approximates the precise numerical 
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FIG. 6. (color online). Off-diagonal elements p„m = {n\pss\m) of 
the steady-states pss as a function of the tuning parameter k for the 
same parameters as in Fig. 2(a). These nonzero coherences demon- 
strate that the steady states are not maximally mixed even when k is 
not an integer. 



evolution of the infinite-dimensional system. This effect we 
refer to as the nonstationary-state single-PB, but it is usually 
called the single-photon optical truncation |4]. 

Finally, we note that these PB effects can also be interpreted 
as photon-induced tunneling studied in the context of PB in, 
e.g., Refs. ifTll [161 [Slll . An especially simple interpretation 
can be found for the dissipation-free PBs. Specifically, we can 
say that the PB effect for k — 2, as shown in Fig. 4, mainly 
corresponds to the dominant transition from the ground to the 
second-excited state with the two-photon resonant transition 
condition as schematically presented in Fig. 1. Analogously, 
the PB effect for fc = 3, as shown in Fig. 5, describes mainly 
the three-photon resonant transition from the ground to the 
third-excited state, which can be explained with the help of 
Fig. 1. Due to amplitude dissipation, the population of lower- 
excited states increases, as seen by comparing Figs. 2(a) and 4 
for the two-PB, as well as Figs. 2(b) and 5 for the three-PB. 



IV. QUANTUM SIGNATURES OF PHOTON BLOCKADES 
A. Photon-number signatures of photon blockades 

We have already discussed some photon-number signa- 
tures of the PB effects by analyzing Fig. 2. Now, we focus 
on demonstrating the nonclassicality of the generated steady 
states. We recall that the nonclassicality (or quantumness) of 
a bosonic system is usually understood if the system is de- 
scribed by a non-positive Glauber-Sudarshan quasiprobability 
function |32] or, equivalently, by a negative normally-ordered 
matrix of moments ll33[|34ll . In particular, the nonclassicality 
can often be revealed by analyzing photon-number properties 
only. 
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Thus, in Fig. 2, we plotted the Fano factor, which is defined 
byil: 



FiPss) = 



(25) 



A given state exhibits sub-Poisson (super-Poisson) photon- 
number statistics if F < 1 (F > 1). The sub-Poisson statistics 
is a nonclassical effect. 

It is seen that the generated state exhibits the sub-Poisson 
photon-number statistics for k < 1.5 and near A: = 2 as 
described by the Fano factors for the parameters chosen in 
Fig. 2(a). Global and local minima of the sub-Poisson statis- 
tics occur for the single-PB (fc = 1) and two-PB (k ~ 2), re- 
spectively. Note the occurrence of local maxima of the Fano 
factors at points relatively close to A; = 2. There are also local 
maxima at fc = 3 and 4 for the parameters e, %, and 7 cho- 
sen in Fig. 2(a). However, for the approximately twice-larger 
driving strength e as in, e.g.. Fig. 2(b), we can also observe 
at k — 3 a local minimum below 1 of the Fano factor, which 
corresponds to the sub-Poisson statistics of the three-PB. 



B. Coherence and entropic signatures of photon blockades 

In order to show how the steady-state solutions pss depend 
on the tuning parameter k, we also analyze their coherence 
properties. 

Off-diagonal elements pnm — (?^|Pss|"^) (for n ^ m) 
of the density matrix p^s, which are also often called coher- 
ences, are shown in Fig. 6. It is seen that almost all plotted 
coherences do not vanish in the steady states even for off- 
resonances (i.e., when the tuning parameter k is not an in- 
teger). This means that the generated steady states are not 
completely mixed. 

Our deeper analysis of the coherence properties of the 
steady states presented in Fig. 7 includes the following mea- 
sures: the purity p{pss) = tr [pl^] and the coherence parame- 
ter C{pss) defined as the sum of all off-diagonal terms of the 
density matrix pss (see, e.g., llssll ): 

C{pss) = ^ IPnmP tr [ipss " Pdiag)^] 

= m(Pss) - Ai(/5diag), (26) 

where pdiag = J2n P™i\n){n\. Thus, this parameter is just 
the total-state purity after subtracting the diagonal-state pu- 
rity. It is worth noting that C is sometimes additionally nor- 
malized [35]. It is seen that decoherence, understood as the 
vanishing of the off-diagonal elements of a density matrix, 
can intuitively be quantified by C{pss)- 

To describe the mixedness of the generated steady state, we 
calculate the von Neumann entropy S{pss) = — tr {pss In Pss) 
and the thermalization parameter defined for a finite- 
dimensional system as ifssll : 

Sl{Pss) 



T{Pss) = 



^/tr {p - /5o)tr (p - pn 
-SlCPs; 
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FIG. 7. (color online). The purity fj,{pss), von Neumann entropy 
S{pss), thermalization Tthcimipss), and coherence C(pss) as a func- 
tion of the tuning parameter k for the mean thermal-photon numbers 
(a) fith = and (b) nth = 0.05 with other parameters chosen the 
same as in Fig. 2(a). This figure shows how thermal photons strongly 
affect p, S, and Tthcrm, but do not affect much C. 



which is defined as 
S'l(Pss)- Here, po = 

(0|pss|0), Pmax = {'' 



a properly normalized linear entropy 

^ |0)(0|, Pinax = |"max)("-max|, PO = 



+ P{Pss) - 2po][l + P(Pss) - 2p 
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= ("^maxIPssl'^max), and lumax) is the Up- 
permost Fock state generated in the system. As explained 
in Ref. fsS*], this normalization of the linear entropy is done 
to exclude somehow the contribution of the ground state 
Po = |0)(0|, which has double nature: It is both a pure 
state and a completely decoherent equilibrium state (as in our 
model without the driving force). Moreover, the linear en- 
tropy Sl{Pss), which is another parameter of mixedness, can 
easily be obtained from the purity (depicted in Fig. 7), since 
Sl{Pss) = 1 — p(Pss)- Other aspects of the quantum entropies 
and mixedness in the discussed model (but only in the special 
(27) case of fc = 1) were studied in Ref. [36]. 

One can clearly see in Fig. 7 that the entropy S{pss) and 
thermalization Tthcrm (Pss) reach maxima and, equivalently, 
' the purity p{pss) and coherence C(pss) parameters have min- 
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FIG. 9. (color online). Spectra of squeezing (a) 50(1^) and (b) 
Stt/2{'-^) for the single-photon blockade for the same parameters as 
in Fig. 2(a). Broken curves show the spectra for the Hamiltonian 
^rot (0) ^''■'^ off-resonance values of the tuning parameter k. 
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FIG. 8. (color online). Wigner functions W{a = x + ip) for the 
fc-photon blockades corresponding to the steady-state solutions of 
the master equation for the Hamiltonian Hj:^l{0), with the same pa- 
rameters as in Fig. 2(a) except the driving strength e = 11. 567 for 
= 3. It is seen that the Wigner function for the fc-photon blockade 
has k peaks and k dips. For A: = 3 the peaks are deformed (due to 
interference in phase space) but still are visible. 



FIG. 10. (color online). Spectra of squeezing for the two-photon 
blockade analogously to those in Fig. 9 but for the resonance fc = 2 
and off-resonance values of the tuning parameter k. It is seen that 
either So or 5.^/2 only for the off-resonant cases has a clear negative 
dip. 



C. Phase-space description of photon-blockades 

The dynamics of quantum systems can equivalently be de- 
scribed by using the phase-space formaUsm of Wigner func- 
tions. This formalism is particularly useful for distinguishing 
different PBs effects as we show below. The Wigner function 
can be defined as |37fl : 



ima at (or very close to) k — 1, 2, 3. It is worth comparing 
this behavior with the Fano factor F{pss) shown in Fig. 2(a) 
for the same parameters, where maxima are observed only for 
k = 1 and 2 (where the PBs are predicted), while for fc = 3 
the minimum occurs, which corresponds to the case where 
the PB is practically not observed for the chosen parameters, 
as we conclude by analyzing the probabilities P„ shown in 
Fig. 2(a). By contrast, the coherence parameter C(/5) has min- 
ima at fc w 1,2 and very small maximum at fc = 3. By com- 
paring Figs. 7(a) and 7(b) obtained for the system coupled to 
a zero- and nonzero-temperature reservoirs, we can easily see 
how the PBs are sensitive to temperature or, equivalently, to 
the number fith of thermal photons. Even by adding a very 
small number of thermal photons, like nth = 0.05, the linear 
and von Neumann entropies together with the thermalization 
parameter are noticeably increased. This is not the case for 
the coherence parameter C(pss)- 



W(a) ^ -tr \b-\a)pD(a)P] , (28) 

TT 

where D{a) = exp(cka^ — a*a) is the displacement operator 
with a complex number a, and P = exp(i7rata) is the par- 
ity operator, so its action on Fock states is simply given by 
P\n) — (— The Wigner function can be generalized 
to the s-parametrized Cahill-Glauber quasiprobabilities [37*|. 
Nevertheless, contrary to other definitions, Eq. (|2FI shows a 
direct method (i.e., without the necessity of applying quan- 
tum state tomography) to measure the Wigner function |38']. 
This direct method was experimentally applied in, e.g., cavity 
QED 1 39] and circuit QED |40] systems. 

In Fig. 8, we show the Wigner function for the blockades 
of up to fc photons for fc — 1, 2, 3 by properly choosing pa- 
rameters in order to satisfy the resonance condition A^: = 
together with 7 ^ e ^ x- The fc-peak and fc-dip (anti-peak) 
structures of the Wigner functions clearly correspond to the 
fc-PB. 
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FIG. 11. (color online). Spectra of squeezing for the three-photon 
blockade analogous to those in Figs. 9-10 but for the resonance fc = 3 
and off-resonance values of the tuning parameter k (assuming the 
driving strength e = 11. 667). Analogously to the two-photon block- 
ade, the off-resonance spectra have a negative dip either for 6 = or 
7r/2, which does not appear for the resonant case for = 3. 



D. Spectrum of squeezing for photon-blockades 

The PBs can also be revealed in two-time correlations. 
Let us analyze the two-time normally-ordered (::) and time- 
ordered (T) correlation function 

r(: Xg{T)XeiO) :) = lim {Xg{t + T)Xe{t)) 

= {a\T)am + {a\Q)a{T)) 

+ e-^'\a[T)a{Q)) + e''\a\G)a) {t)) (29) 

of the quadrature phase distribution 

Xg{t) = a{t)e-'^ + a){t)e'^. (30) 

The spectrum of squeezing is defined as the Fourier transform 
of the covariance 114111 : 

/oo 
dTe~'^^T{:Xe{T),Xe{Q):), (31) 
-00 

where the covariance is defined by the general formula 
{A,B) = (AB) - (A) (B). In Figs. 9-11, we demonstrated 
distinctive properties of the spectra of squeezing for the 
single-PB (Fig. 9), two-PB (Fig. 10) and three-PB (Fig. 11) 
for the tuning parameters at the resonance, k = fco. and 
slightly out of the resonance, k = ko ± 6, where (5 <C 1 and 
ko ~ 1, 2, 3. Surprisingly, the spectra of squeezing as a func- 
tion of the tuning parameter k are more sensitive indicators of 
the two- and three-PBs instead of those of the single-PB. 

V. CONCLUSIONS 

We studied nonlinear photon-photon interaction at the 
single-photon level in a system consisting of a cavity with 
a Kerr nonlinearity driven by a weak classical field. This is 
a standard prototype model, where the KeiT-like nonlineaiity 
can be induced by the interaction of, e.g., a qubit with a cavity 
field in the dispersive regime. 

By finding the master-equation solutions of the model in the 
steady-state limit for a properly chosen frequency of the clas- 
sical driving field, we observed a blockade of more than one 



photon. This is a generalization of the single-PB for the multi- 
photon (say TO-photon) case, which means that the Fock states 
with n = 0, 1, m photons are, practically, the only gener- 
ated in this nonlinear cavity. This effect can also be inter- 
preted as the multiphoton-state truncation realized by nonlin- 
ear quantum scissors or multiphoton-induced tunneling corre- 
sponding to multi-photon transitions as schematically shown 
in Fig. 1. 

It is worth noting that some generalizations of the single- 
photon nonlinear tiTincation processes for the multiphoton 
cases were already discussed by, e.g., Leoriski et al. (see 
Ref. |4] for a review). Nevertheless, all these proposals as- 
sume either higher-order driving processes, as described in 
ReL iH by the Hamiltonian e(')[a' + {a)% or 

higher-order Kerr nonlinearity, corresponding to the Hamilto- 
nian hI^L = x^'"^^Ha^)''^^a''+^ in Ref. |23], or both ef'^) 
andx^'^ II22II ioikj = 2,3,.... In our approach, we assume 
the lowest-order parametric driving (/ = 1) and the lowest- 
order Kerr nonlineaiity (fc ~ 1), which is the same as in the 
standard single-PB. Also, in the mentioned generalizations, 
only the free evolution was analyzed. Thus, these effects can 
solely be interpreted as nonstationary-field effects; so, they are 
almost unmeasurable. The dissipation effects were later stud- 
ied, but only for the short-time evolution regime in Ref. [:4i. 
In this paper, we focused on measurable effects, thus, we stud- 
ied (except only one section) the steady-state solutions of the 
master equation in the infinite-time limit. 

We described, as given by Eq. (|9]l, how to choose the reso- 
nance frequency Ud of the driving field for a given cavity fre- 
quency ujQ and the KeiT nonlinearity x and how to increase the 
driving strength e (see Fig. 3) in order to observe the blockade 
at the multiphoton Fock state. 

In our presentation we focused on the blockades of fc = 
2 and 3 photons in order to achieve high fidelity Fk of the 
fc-photon truncation and with a contribution of the fc-photon 
state which cannot be neglected, which means that ^ 1. 

We showed analytically and numerically that these PBs oc- 
cur both for non-dissipative and dissipative systems. In par- 
ticular, we showed an excellent agreement (i.e., no differences 
on the scale of Fig. 4) for the two-PB between the numerical 
solutions obtained in the Hilbert spaces of infinite-dimension 
(practically A^dim — 100) and of finite dimension A^dim = 4). 
Our approximate analytical solutions only slightly differ from 
the numerical ones (as shown in Fig. 4 for a relatively large 6). 

We demonstrated a vaiiety of quantum properties reveal- 
ing the unique nature of the single-, two- and three-PBs as a 
function of the tuning parameter fc. In particular, we stud- 
ied photon-number statistics (as shown in Fig. 2), as well 
as coherence and entropic properties (in Figs. 6 and 7). We 
gave a clear comparison of the Wigner functions of the steady 
states corresponding to the blockades of fc = 1,2,3 photons 
(see Fig. 8). Moreover, we showed clear differences in the 
spectra of squeezing of the steady states in the resonant cases 
at fc = 1,2,3 and slightly off the resonances (as shown in 
Figs. 9-11). 

We suggested that the two-PB and three-PB can be ob- 
served in various systems, where the single-PB has already 
been experimentally observed lfT3l[l4l[T6ll or, at least, theoret- 
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ically predicted flT^ . 

Analogously, a blockade of phonons instead of photons 
can be considered. Namely, a two-phonon generalization of 
the standard single-phonon blockade can be predicted in the 
nanomechanical systems studied in Refs. [18,, 191 . The cru- 
cial point is that the effective Hamiltonians, given by Eqs. ([T]i 
and under proper conditions, can be obtained from other 
standard models in cavity QED, circuit QED, and quantum 
optomechanics. 

The single-photon blockade has attracted considerable in- 
terest, with potential applications in quantum state engineer- 
ing, quantum information, and quantum communication. We 
believe that the multi -photon blockades described in this paper 
can also find some useful quantum applications and also, at a 
fundamental level, they can show a deeper analogy between 
condensed-matter and optical phenomena. 
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